Butterfly Floquet Spectrum in Driven SU(2) Systems 
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The Floquet spectrum of a class of driven SU(2) systems is shown to display a butterfly pattern 
with multi-fractal properties. The level crossing between Floquet states of the same parity or dif- 
ferent parities is studied. The results are relevant to studies of fractal statistics, quantum chaos, 
coherent destruction of tunneling, and the validity of mean-field descriptions of Bose-Einstein con- 
densates. 
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Hofstadter's butterfly spectrum of the Harper model 
[l[ has attracted tremendous mathematical, theoretical 
and experimental interests. For an arbitrary irrational 
value of one system parameter, the spectrum of the 
Harper model is a fractal, which has been strictly proved 
after decades of research on the "Ten Martini problem" 
As one important implication, a fractal butterfly 
spectrum suggests the closing of a quantum gap inflnite 
times and hence the occurrence of infinite quantum phase 
transitions 3]. 

Early quantum chaos studies established that the Flo- 
quet (quasi-energy) spectrum of periodically driven sys- 
tems may display a fractal butterfly pattern as well 
However, the nature of the fractal Floquet spectrum is 
still poorly understood. Indeed, because the eigen-phase 
of Floquet states is restricted to a range of 2tt, under- 
standing a Floquet spectrum associated with an infinite- 
dimensional Hilbert space is subtle and challenging 
Furthermore, previous findings regarding to the fractal 
Floquet spectrum were largely limited to the so-called 
kicked-Harper model (a driven version of the Harper 
model) [7, 8, 9] and its variant [HI. [ol. \ld^ . 
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promising possibilities for the engineering and simula- 
tion of driven ultracold systems with a prescribed Flo- 
quet spectrum. In this Letter, we show that the Flo- 
quet spectrum of a deceptively simple class of SU(2) sys- 
tems, constructed from a driven two-mode Bose-Einstein 
condensate (BEG), displays a butterfly pattern and pos- 
sesses truly remarkable properties. For example, we show 
that with one certain system parameter fixed the overall 
butterfly pattern is insensitive to the number of bosons 
(denoted N) in the BEG, but some detailed features de- 
pend on whether N is odd or even. We shall reveal that 
the found butterfly pattern contains many level cross- 
ings between states of different parities and thus many 
points of coherent destruction of tunneling (GDT) [15j . 
with the total number of GDT points found to scale as 



~ N^ '^. As an analog of first-order quantum phase tran- 
sitions, we discover that the found butterfiy pattern also 
contains many level crossings between same-parity eigen- 
states. These results suggest that the class of driven 
SU(2) systems studied here may become a test bed for a 
number of research topics. Several specific applications 
of this work are also discussed. 

Driven two-mode BEG systems were proposed before 
11 [H [ill to realize the well-known kicked top model 
18| in the quantum chaos literature. In its most general 
form, a driven two-mode Bose-Hubbard model can be 
written as 



H = f{t)h{ala2 



ajfli) 



+ g{t)h,{a\ai - ala2)^, (1) 



where and a- are the bosonic annihilation and cre- 
ation operators for the ith mode, f{t) describes the time- 
dependent tunneling rate between the two modes, and 
the g{t) term describes the self-interaction between same- 
site bosons, whose time-dependence can be achieved by 
Feshbach resonance induced by a magnetic field. The 
total number of bosons N — a\ai + a\a2 is a con- 
served quantity and the dimension of the Hilbert space 
is A'^ -|- 1. Using the Schwinger representation of angu- 
lar momentum operators, namely, Jx — (a|a2 -t'a2ai)/2, 
J„ = (ajfli — a|a2)/(2z), and — (ajai — 0302)72, Eq. 



dl]) reduces to 



H = 2fit)nJx + 4,g{t)hJ^ 



(2) 



Glearly, the dynamics is solely determined by the SU(2) 
generators J^, Jy and Jz- The total angular momentum 
quantum number J is given by J = N/2. The Hilbert 
space can be expanded by the eigenstates of Jz, denoted 
|m), with Jz|m) = m\m). The population difference be- 
tween the two modes is given by the expectation value of 
2Jz. If we exchange the indices of the two modes, then Jx 
is invariant, Jz —>■ — Jz, and as a result the Hamiltonian 
in Eq. ^ is unchanged. This refiects a parity symmetry 
of our model. 
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Consider then two specific forms of f{t) and g{t). In 
the first case f{t) = a/{2T), g[t) = goY^J[6{t - 2nT - 
t) — 5{t — 2nr)]. The Floquet operator, i.e., the unitary 
evolution operator F from 2nr + 0+ to (2?! + 2)t + 0+, 
is then given by 

where ry = 4:goN. Because the first two or the last 
two factors in Eq. ([3]) constitute the Floquet oper- 
ator for a standard kicked-top model [l8], our driven 
system here can be regarded as a "double-kicked top 
model". Alternatively, if we set g{t) = go/S,, J{t) = 
f +'^(^ ~ where ^ is the time delay 

between the two delta kicking sequences, then the asso- 
ciated propagator F' from nr — 0"*" to {n + \)t ~ 0+ is 
given by 

pi ^ g4J?jf/2,/g-i(4gor/C)jJg-4aJxg-i»)Jj/2Jg-iQJx_ (-4-) 

Under the special condition 4(7ot/^ = 2fc7r (S/ctt) for inte- 
ger J (half integer J), where k is an integer , the factor 
g-«(4go-r/C)J3 jg yjij-ty in the (2 J -I- l)-dimensional Hilbert 
space and hence F' becomes identical with F. Thus, 
there exist two different scenarios for realizing F, the 
key operator to be analyzed below. 

In the |m) representation, the third factor e""''^^/^'^ 
of F equals e^*''™ which is a pseudo-random num- 
ber for irrational "q/J. The first factor of F however 
effectively induces a time-reversal of the third factor and 
thus partially cancels this pseudo-random phase. Indeed, 
using the SU(2) algebra, the product of the first three 
factors of F in Eq. ([3]) is given by 

^ g-»c{(J./2+.J„/2)e*I''(=^-+i)/=''l+c.c.}^ (5-) 

showing that the 77-dependent term entering into F be- 
comes e*t'''^^'^^+^)/^'^l , which is always a quasi-periodic 
number ^'■['n{'i™.+i)/2j\ ^^^q \rn) representation. This 
partial cancelation of quasi-random dynamical phases 
implies intriguing spectral properties [19| . 

To study the classical limit of F we consider scaled 
variables x = Jxl J , y = Jy/J, and z = Jz/J- The three 
operators a:, y, z also satisfy the angular-momentum al- 
gebra, but with an effective Planck constant ?ioft = 1/J. 
Taking the ?icff — > limit with fixed 77 and a, the clas- 
sical dynamics associated with F can be obtained, with 
variables x, y, and z restricted on a unit sphere. Be- 
cause 77 = 4,goN, this classical limit with fixed 77 requires 
N — > -|-oo and go ^ 0. This condition is apparently 
equivalent to that in a standard mean-field limit of the 
driven BEC. 

Figure 1 shows the typical eigen-phase spectrum of F 
vs = rjhcff = 7]/ J = 8go, for J = 20,20.5,100 and 
a/hes — 1.0. Because the spectrum of F is invariant 
if fir^ hjj + Att, we set hjj G [0,Att). Though in Fig. 




FIG. 1: Eigen-phase spectrum (denoted e) of the Floquet 
operator F in Eq. (3). J = 20 in (a), 20.5 in (b), and 100 in 
(c). a/hcB — a ■ J = 1.0 in all panels. 

1 the involved Hilbert space is rather small, spectacular 
butterfly patterns are already obtained (their symmetry 
with respect to — 2t: can be proved). They resemble 
the famous Hofstadter's butterfly, but also present re- 
markable differences in several aspects. First, if we take 
a vertical cut of the butterfly patterns in Fig. 1, the spec- 
trum is not found to present any large gap. Second, the 
butterfly patterns shown in each panel of Fig. 1 possess 
a double-butterfly structure, with each butterfly cover- 
ing a 27r range of This double-butterfly structure 
is somewhat analogous to the spectrum of a Harper-like 
effective Hamiltonian considered in Ref. [Q]. More in- 
terestingly, though Fig. 1(c) has much more levels than 
Fig. 1(a) and Fig. 1(b), the overall outline of the double- 
butterfly structure is seen to be insensitive to J for fixed 
Q:/?icff — a ■ J. For a fixed value of J but for other not 
too large values of a, the qualitative features of the but- 
terfiy spectrum remain, but at different scales. For very 
large values of a (e.g., a/hcs > 10), the butterfly pattern 
for a fixed value of J will gradually dissolve, in a similar 
fashion as in the kicked Harper model 4;]. 

Some detailed aspects of the spectrum are also note- 
worthy. For example, it is observed that the spectrum 
collapses to one point for hrj = 27r, if and only if J is 
an integer. This can be explained as follows. If J is an 
integer and if = 27r, then in the |m) representation, 

^-ivJl/2J ^ g-iTrm^ ^ g-iTrm ^ 5"""''=. So in this CaSC 

Q-^ri'h/'^J is equivalent to a rotation of tt around the z 
axis, and hence the first three factors of F exactly cancel 
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FIG. 2: Poincare surfaces of section (with J^, > 0) of the 
classical or mean-field limit of F in Eq. (3), with q = 0.05 
(same as in Fig. 1(a)), r; = 5 in (a), 30 in (b), and 75 in (c). 



its last factor. This cancelation will not occur if J is a 
half integer, i.e., if N is odd. 

Figure 2 depicts the phase space structure of the clas- 
sical limit of F. As 77 increases, the classical dynam- 
ics changes from being regular to being chaotic. On the 
other hand, the Floquet spectrum shown in Fig. 1 can 
be however much similar for radically different values of 
77 ranging from to 47r/?ioff — 4j7r. Therefore, upon 
quantization the regular or chaotic nature of the clas- 
sical dynamics might not necessarily be reflected in the 
spectrum and hence can be irrelevant to the quantum 
dynamics. 

The statistical behavior of the found butterfly spec- 
trum is also examined. To have good statistics we con- 
sider a much larger value of J. Figure 3(a) presents the 
cumulative level density TSS(e) for a representative value 
of hrf. N{e) is highly irregular, but does not show any 
clear flat steps. This is consistent with our early obser- 
vation that no large gap exists in the spectrum. Figure 
3(b)-3(d) shows the associated level distribution P(e) at 
three different scales. Evidently, P{e) has a fascinating 
self-similar property. This motivates us to quantitatively 
characterize the spectrum via the generalized fractal di- 
mension Dq, with the results shown in Fig. 3(e). As 
expected from the N{e) result in Fig. 3(a), Do = 1. 
However Dq for q ^ clearly shows that the spectrum 
has multi-fractal properties. For comparison, Fig. 3(e) 
also shows the Dq result for a standard kicked-top model 
with the same values of rj and a (i.e., considering an op- 
erator comprising only the first two factors of F). The 
Dq behavior in the kicked-top case is as trivial as that 
of a random sequence: it remains close to unity and 
slightly decreases with increasing q due to finite-size ef- 
fects. Based on these results, we conjecture and invite 
a formal mathematical proof that the butterfly patterns 
found here contain true fractals in the limit of J -l-oo. 

We next study the level crossings in the butterfly pat- 
terns. Interestingly, the minimal distance in between 
two level crossings is found to decrease sharply with J. 
So even for a rather small J 10 it is already computa- 
tionally demanding to identify all the level crossings. As 
an example Fig. 4(a) presents the typical level crossing 
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FIG. 3; The cumulative Floquet state density (a) and the 
Floquet state density distribution (b)-(d) at different scales, 
for hr, = {V5 - l)7r/2, a/h^a = 1, and J = 2999. Panel 
(e) shows the generalized fractal dimension Dq. Crosses and 
circles are for odd-parity and even-parity states. Triangles 
represent the result for a standard kicked-top model. 



behavior in the vicinity of a null cigcn-phase for J — 10. 
The Floquet states are seen to cross each other frequently, 
between different-parity states and between same-parity 
states. Both types of level-crossings are of enormous 
interest. For the first type, at a crossing point an ar- 
bitrary superposition of two crossing states of different 
parities remains an eigenstate but generally breaks the 
parity symmetry, thereby maintaining a nonzero popula- 
tion difference between the two modes forever [13|. This 
makes it clear that the first type of level crossings give rise 
to the seminal CDT phenomenon 15| that has attracted 
broad experimental interests. Note that in some regimes 



of h,^ , to the naked eye two curves of opposite parities in 
Fig. 4(a) are almost on top of each other, and as a result 
many CDT points are found in these regimes. Note also 
that CDT-induced population trapping is fundamentally 
different from the well-known self-trapping effect on the 
mean-field level. Indeed, the CDT effect here depends on 
rj and J, whereas mean-field self-trapping is transient and 
independent of J. Now turning to the second type of level 
crossings, they come as a surprise because avoided cross- 
ings between same-parity states are generally anticipated 
for classically non-integrable systems (see Fig. 2). The 
second type of crossings therefore suggest the uniqueness 
(e.g., some effective local "symmetry") of F whose ma- 
trix elements in the |m) representation are quasi-periodic. 
Recalling the above-mentioned extreme example where 
all levels cross at hr, = 2it for integer J, we expect that 
special arithmetic properties of play a key role in both 
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FIG. 4: (a) Level crossings between three even-parity states 
(dashed lines) and two odd-parity states (solid lines), for J = 
10 and a/hefi = 1.0. (b) Number of level crossings versus J, 
for £ [0, 47r) and a/heS = 1-0. The cross (square) symbols 
are for crossings between different-parity (same-parity) states 
and the fitting suggests a power law scaling J^" (J^'^). 



types of level crossings. 

By obtaining all the level crossings with high accuracy 
for J < 12, we obtain in Fig. 4(b) that the number of 
CDT points contained in the butterfly patterns scales as 
J'^" and the number of same-parity crossings scales as 
J^-^. In particular, we conclude that as N goes to infin- 
ity, on average each pair of Floquet states in a butterfly 
pattern see inflnite CDT points. 

In the kicked Harper model, the quantization rule 
varies with the boundary condition adopted [20] and a 
compact toroidal phase space arises only if the Planck 
constant assumes special values [ilj. A general treat- 
ment of the kicked Harper model leads to a band struc- 
ture that often complicates the issue. By contrast, the 
phase space here is necessarily on a sphere [l^, with 
no arbitrariness in quantization and no band structure 
in the spectrum. For these reasons the new butterfly 
Floquet spectrum discovered in this work can stimulate 
more studies of fractal Floquet spectrum in driven sys- 
tems. Results here also suggest that our strategy in gen- 
erating a butterfly quasi-energy spectrum, namely, the 
use of partial cancelation of quasi-random phases (first 
advocated in a kicked-rotor system [5]), is widely appli- 
cable. Furthermore, it is now clear that three quantum 
chaos paradigms, i.e., the kicked-top, kicked-Harper, and 
kicked-rotor models, are all linked together for the first 
time, because all of them can display fractal statistics. 

Finally we mention two specific applications. First, be- 
cause the found butterfly spectrum collapses at Tiri = 27r 
(or go = 7''/4) for integer J, one may experimentally de- 
termine if N is even or odd by scanning the dynamics 
in the neighborhood of go = it /A. This possibility does 
not exist in the mean-field dynamics of a BEC. Similarly, 
one may study the CDT points to reveal non-mean-field 



effects. Second, it is now of great interest, both experi- 
mentally and computationally, to revisit early results of 
how a multi-fractal spectrum can be manifested in time- 
dependent properties 

Detailed results of this work will be published else- 
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by WBS grant Nos. R-144-050-193-101/133 and R-144- 
000-195-123. 
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